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Plasmons produce large confinement and enhancement of light that enable applications as varied 
as cancer therapy and catalysis. Adding to these appealing properties, graphene has emerged as 
a robust, electrically tunable material exhibiting plasmons that strongly depend on the density of 
doping charges. Here we show that adding a single electron to a graphene nanoisland consisting of 
hundreds or thousands of atoms switches on infrared plasmons that were previously absent from the 
uncharged structure. Remarkably, the addition of each further electron produces a dramatic fre- 
quency shift. Plasmons in these islands are shown to be tunable down to near infrared wavelengths. 
These phenomena are highly sensitive to carbon edges. Specifically, armchair nanotriangles display 
sharp plasmons that are associated with intense near- field enhancement, as well as absorption cross- 
sections exceeding the geometrical area occupied by the graphene. In contrast, zigzag triangles do 
not support these plasmons. Our conclusions rely on realistic quantum-mechanical calculations, 
which are in ostensible disagreement with classical electromagnetic simulations, thus revealing the 
quantum nature of the plasmons. This study shows a high sensitivity of graphene nanoislands 
to elementary charges, therefore emphasizing their great potential for novel nano-optoelectronics 
applications. 



I. INTRODUCTION 

Light can efficiently excite plasmons (collective elec- 
tron oscillations in matter), thus producing strong con- 
finement of electromagnetic energy and huge enhance- 
ment of the associated electric fields. These phenomena 
have spurred a plethora of applications ranging from ul- 
trasensitive detection p~H5] to cancer therapy [6j [7] and 
catalysis [8 , 9 , which have been largely fueled by progress 
in the synthesis of noble-metal nanoparticles with in- 
creasing control over size and morphology [10] . Plasmons 
in conventional metals result from the cooperative effect 
of many electron-hole (e-h) virtual excitations around the 
Fermi level. The resulting plasmon frequencies scale with 
the 1/2 power of the density of valence electrons. Under- 
standably, massive amounts of charging are required to 
produce sizable frequency shifts in the plasmons (TTJ [12] , 
and therefore, their electrical control remains elusive. 

This scenario has substantially changed with the ar- 
rival of graphene. Plasmons exist in this material when 
it is doped, but now the plasmon frequency is propor- 
tional to the 1/4 power of the doping charge density 
rather than the square root of the valence electron den- 
sity [13]. This behavior is a consequence of the peculiar 
electronic structure of graphene, characterized by a van- 
ishing density of electron states at the Fermi level. Evi- 
dence for graphene plasmons has been recently reported 
through terahertz [14] and infrared (IR) [T5UT7] optical 
spectroscopies, as well as through direct near-field spa- 
tial imaging [18] [19]. These studies have conclusively 
demonstrated that graphene plasmons can be frequency- 
tuned using conventional electric gating technology. In a 
parallel promising effort, the electro-optical tunability of 
graphene has also been used to modulate the plasmonic 



* Electronic address: J.G.deAbajo@nanophotonics.es 



response of neighboring metallic nanoparticles [2QU22] . 
Now, the question arises, is the singular electronic struc- 
ture of extended graphene permeating the optical prop- 
erties of nanometer-sized doped graphene islands as well? 
Can we bring the level of doping needed to sustain plas- 
mons in a small structure down to just a single electron? 

Here we predict that IR and near-IR (NIR) plasmons 
in graphene nanoislands can be switched on and off by 
the addition or removal of a single electron. This is a 
remarkable property considering that the islands contain 
hundreds of atoms. Specifically, we focus on graphene 
nanotriangles, which we describe using a combination of 
a tight-binding model for the electronic structure and the 
random-phase approximation (RPA) [23] for the dielec- 
tric response, as explained in Appendix [A] A more de- 
tailed description of the specifics of our approach is given 
elsewhere [24], essentially extending to finite graphene 
islands what was previously reported for more extended 
systems by combining tight-binding and RPA. This effort 
was pioneered by Wallace [25] in graphene and graphite, 
and continued through outstanding analyses by other au- 
thors [H (21 123 . 



II. RESULTS AND DISCUSSION 

Figure [I] summarizes our main finding: an electrically 
neutral graphene nanotriangle (with armchair edges and 
7nm side length) does not display strong polarization 
when illuminated by IR light; in contrast, the addition 
of a single electron results in the emergence of an intense 
plasmon mode, so that strong polarization is produced 
when the light is tuned to the plasmon energy (0.38 eV). 
A single electron can thus trigger the existence of plas- 
mons in the structure. When further electrons are added 
to such armchair nanotriangle (Fig. [2^), this IR plasmon 
undergoes a dramatic blue shift (in steps of ~ 0.1 eV) 
and it gains in strength. Notice that we denote the net 
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FIG. 1: Single-electron switching of graphene plasmons. A neutral graphene nanotriangle (left) shows negligible polar- 
ization under external illumination. The same nanoisland displays a 0.38 eV plasmon resonance when it is charged with one 
electron, thus undergoing strong polarization (right). The density plots show the induced electron charge at the carbon sites. 
For comparison, the extinction cross-section is a sizable fraction of the graphene area in the resonant charged configuration 
and two order of magnitude smaller in the neutral island. 



charge of the structure as Q in Fig. |2j so that Q < 
(Q > 0) corresponds to doping with electrons (holes). 

Equally important, we find a NIR plasmon with and 
without doping at an energy around ~ 0.88 eV. In con- 
trast to the IR plasmon, this high-energy mode undergoes 
redshift with increasing doping. Both the IR and the 
NIR plasmons give rise to large absorption cross-sections 
(Fig.^), which can reach values exceeding the geomet- 
rical area of the graphene for large doping in the IR and 
for low doping in the NIR (see Appendix [i]) . 

The plasmon shift is substantially larger than the 
width, thus making it clearly resolvable. Obviously, this 
conclusion depends on the parameter u sed f or the intrin- 
sic width Tit -1 in the RPA (see Fig. A2 in Appendix 
0. The main contributions to the width in extended 
graphene originate in optical-phonon losses, impurities, 
and disorder. These mechanisms are well described in 
extended graphene [28], and we take fir -1 = 1.6 meV 
(i.e., a dephasing time r ~ 400 fs) as a reasonable es- 
timate for high-quality samples [29, 30]. Actually, given 
the small area of the islands, it should be possible to iden- 
tify many of them without impurities or defects. Besides, 
the large shift-to-width ratios of Fig. [2] guarantee that 
our conclusions are still maintained with much higher 
losses up to the upper limit that is intuitively imposed 
by the lifetime of hot electrons (~ 85 fs, as resolved from 
two-photon photoemission [31]). Additionally, Landau- 
damping is expected to be negligible for low-energy plas- 
mons in defect-free extended graphene [32 ] , and although 
plasmon decay into e-h pairs is made possible by the loss 
of translational symmetry in nanoislands, no substantial 
broadening is observed in the plasmons of armchair trian- 
gles beyond the intrinsic width Tit -1 introduced through 



the RPA formula (see Fig. |A2| and discussion below). Ac- 
tually, the plasmon energies do not overlap any intense 
e-h transition (see Fig. [3^ below). 

Incidentally, a classical electromagnetic description of 
the graphene (see Appendix |A| also yields IR plasmons 
with similar shifts as the RPA calculations (Fig. [2|l, blue 
curve), although the plasmon energies are substantially 
higher in the latter due to quantum confinement [24] . 
However, the NIR plasmons and their redshifts with in- 
creasing doping are completely missed by classical theory. 

A completely different scenario is observed in zigzag 
nanotriangles (Fig. [2j),d) compared to armchair struc- 
tures (Fig. [2^l, c): a single plasmon appears at ~ 0.78 eV 
instead of the IR and NIR plasmons; and the plasmon 
energy is rather independent of doping. This seems to be 
connected to the presence of near-zero-energy electronic 
states, which are known to exist near zigzag edges [33]- 
[38] . In a previous study [24], these states were found to 
produce plasmon dephasing in small islands. In the nan- 
otriangle of Fig. J2|d, the number of such states (24 per 
electron spin, see Appendix[F| is large enough so that the 
addition of extra electrons does not substantially change 
the electronic structure. Moreover, these states do not 
produce dipole transitions to higher-energy states (see 
Fig. |3]i), and therefore, they play a dummy role in the 
formation of plasmons. For these reasons, mid-IR plas- 
mons are not observed in the zigzag island, and the NIR 
plasmon is rather insensitive to the number of electrons 
added to it, until the zero-energy electronic levels are 
completely filled (plasmons shifts are predicted beyond 
this level of doping, as shown in Appendix 

In extended graphene, the Dirac-cone band structure 
leads to a gap between intra- and inter-band e-h transi- 
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FIG. 2: Plasmon resonances driven by single- and few-electrons or holes in graphene triangular nanoislands. 

We consider nanotriangles of either armchair (ac, a, consisting of N — 816 atoms) or zigzag (zz, b, N = 726) edges (side 
length a = 7nm). c, Absorption spectra of the ac island for various charge states Q/e, obtained from a quantum- mechanical 
description. Doping with electrons (holes) corresponds to Q < (Q > 0). The absorption cross-section (color scale) is 
normalized to the graphene area. The plasmon energies obtained from a classical electromagnetic description of the graphene 
are shown by blue curves, d, Same as c for the zz island, e-h, Near-field intensities (e,f) and charge densities (g,h) associated 
with selected plasmons of c,d for specific values of the doping charge Q/e and the plasmon energy E p . 



tions. Long-lived plasmons exist in that gap. In nanois- 
lands, the electron parallel momentum is not a good 
quantum number due to the lack of translational sym- 
metry, but we can obtain insight into the role played by 
the electronic structure by analyzing the spatial Fourier 



transform of the electron wave functions (Fig. |3pi,b) . An 
incipient Dirac-cone structure is observed despite the fi- 
nite size of the islands, which extends to a well defined 
Dirac cone in larger structures (see Appendix |rl|). In neu- 
tral nanostructures, all states below zero energy are oc- 
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FIG. 3: Understanding plasmons from the electronic structure. a,b, Intensity of the spatial Fourier transform of the 
electron wave functions in the same nanoislands as in Fig. [2] The intensity is summed over all one-electron states and it is 
represented as a function of energy and wave vector along the TK direction. The Dirac cone of extended graphene is shown 
by dashed lines. c,d, Plasmons and e-h excitations mismatch, represented through the single-electron dipole-transition matrix 
elements | {^i\x\^j) | 2 + \(ipi\y\ipj)\ 2 as a function of initial (Ei) and final (Ej) electron-state energies. The area of the symbols 
is proportional to the squared matrix elements. The dark blue lines are obtained by setting \Ei — Ej \ to the plasmon energies 
of the islands under consideration (see Fig. |2t,d). 



cupied, and remarkably, the density of states vanishes at 
the Dirac point in armchair islands (this is because both 
carbon sublattices have the same number of atoms [33 J, 
which explains why the addition of a few electrons causes 
such dramatic changes in the optical response. In con- 
trast, zigzag islands display an intense zero-energy fea- 
ture associated with edge states [33] (see Appendix |H|) , 
so that extra electrons do not produce significant effects 
in the electronic density of states, and therefore, the op- 
tical response is rather insensitive to the net charge of 
the structure until all zero-energy states are occupied. 

The switching on of an IR plasmon with the addition of 
one electron poses the question, is this plasmon consist- 
ing of the oscillation of a single electron? The e-h exci- 
tation spectrum in the neutral armchair island has a gap 
~ 0.7eV, and the plasmon energies are actually not over- 
lapping with those of intense e-h transitions, as we show 



in Fig. [3]3,d. The observed plasmons are thus involving 
multiple interactions among virtual e-h excitations, lead- 
ing to collective electron motion. Further evidence for 
this is obtained by realizing that the observed plasmon 
energies and their characteristics are completely missed 
within a non-interacting electron-gas picture (see Ap- 
pendix [Dj. Nonetheless, despite this mounting evidence 
of multiple e-h interactions, the role of self-screening of 
each electron is still a pending issue that deserves a more 
detailed analysis outside the scope of the present work. 

Plasmon excitations in these graphene islands produce 
a concentration of the electromagnetic energy down to 
a region ~ 10 3 smaller than the wavelength. This en- 
ergy concentration leads to large levels of the field en- 
hancement (~ 10 3 in intensity) upon external illumina- 
tion (Fig.[2^,f), as well as absorption cross-sections larger 
than the geometrical area of the graphene (see Appendix 
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[]| . The induced charge densities associated with these ex- 
citations are clearly showing dipolar excitation patterns 
with finer sign oscillations differentiating plasmons at dif- 
ferent energies (see Fig. |2^,h). 

III. CONCLUSION AND PERSPECTIVES 

We find it remarkable that the addition or removal 
of a single valence electron can switch on and off sharp 
plasmons in a structure already containing hundreds 
of valence electrons. A qualitatively similar conclusion 
is drawn by just examining the response of graphene 
nanoislands from a classical electrodynamics viewpoint. 
However, our quantum-mechanical simulations show that 
this intuition is strongly amended by the details of the 
atomic structure, to the point that the low-energy plas- 
mons under consideration are simply absent from nano- 
triangles with zigzag terminations. In contrast, sharp IR 
plasmons appear in armchair-edge islands, although they 
are severely shifted with respect to classical theory. The 
predicted plasmons, and even their mere existence, are 
thus non-trivially depending on the quantum mechanical 
properties of the underlying graphene fabric. 

An important aspect of electrically tunable 
nanographene plasmons is that they can reach the 
NIR regime due to the reduced size of the structures. 
This is already clear in the armchair triangle of Fig. 
[2^i, but even higher plasmon energies reaching into the 
visible are obtained by reducing the size of the structure 
(see Appendix |C|) . 

The plasmon width in armchair nanotriangles is essen- 
tially limited by the intrinsic decoherence time. Using 
realistic values for this parameter, we predict absorp- 
tion cross-sections exceeding the geometrical area of the 
graphene (see Fig. [2^ and Appendix[C|. This should allow 
the design of patterned graphene sheets with spacings of 
only a few nanometers for complete optical absorption at 
electrically tunable IR and NIR frequencies [39] . 

Graphene nanoislands of sizes in the range of those 
considered here have been fabricated with different meth- 
ods [40-42J. However, a major challenge is the synthe- 
sis of structures with the desired edge terminations. Al- 
though armchair edges are more energetically favorable 
and they are observed in experiments [43j [44], zigzag 
edges grow faster and dominate uncontrolled growth 
[45] [46]. Here, we are predicting appealing optical prop- 
erties for armchair nanoislands, which are nonetheless 
expected to grow stably under very low or very high hy- 
drogen chemical potential [47]. In this regard, a chemical- 
synthesis bottom-up approach can be beneficial to pro- 
duce nanometer-sized graphene structures with well con- 
trolled edges @8j [49] . 

Our prediction of single-electron-induced extreme plas- 
mon shifts relative to the plasmon widths opens new pos- 
sibilities for ultrasmall sensors based upon the observa- 
tion of minute amounts of charge transfer. For exam- 
ple, single-molecule detection should be possible by mea- 



suring the energy shift associated with the transferred 
charge upon absorption of donor or acceptor molecules, 
thus suggesting an alternative optical approach to electri- 
cal single- molecule detection in graphene [50]. However, 
giving the large mismatch between the wavelengths of 
graphene plasmons and light, plasmon readout of indi- 
vidual nanoislands is currently a challenge. In this re- 
spect, electrical [51] or electron-beam plasmon excita- 
tion and detection are promising options. In particu- 
lar, electron beams can maximally couple to the plas- 
mons of graphene. The large wavelength mismatch also 
leads to unprecedentedly high values of the Purcell fac- 
tor (quality factor divided by plasmon mode volume in 
units of the cubed wavelength), the field enhancement 
(with potential application to surface-enhanced IR ab- 
sorption [52 ), and the local density of optical states 
(reaching the strong light-matter interaction regime [53 ). 
Tightly bound plasmons in graphene nanoislands are thus 
a promising tool for the investigation of fundamental op- 
tical phenomena and for applications to sensing and opto- 
electronics. 
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Appendix A: Calculation methods 

1. Quantum- mechanical plasmon calculations 

The linear optical response of graphene nanoislands is 
described in the RPA [23] , using tight-binding electron 
wave functions as input [25]. For the low photon ener- 
gies under consideration, only tt valence electrons con- 
tribute to the response. The electron wave functions are 
expressed in a basis set consisting of one 2p orbital per 
carbon site (denoted |Z)), oriented perpendicularly with 
respect to the graphene plane and occupied with one elec- 
tron on average in undoped structures. A tight-binding 
Hamiltonian is formulated in such basis set, with nonzero 
elements connecting nearest neighbors through a — 2.8eV 
hopping energy [13]. One-electron states and 
energies Sj are obtained upon diagonalization of this 
Hamiltonian, yielding the well-known Dirac-cone elec- 
tronic structure of graphene [13 , 54 . Here, aji is the am- 
plitude of orbital |Z) in state j. For the graphene islands 
and charge states under consideration, electron states are 
filled up to affordable values of the Fermi energy Ep < eV 
relative to the Dirac point. The site-dependent induced 
charge is then obtained from the self-consistent potential 
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is the non-interacting RPA susceptibility, uj is the light 
frequency, fj = {exp[(^£j — Ep)/kBT) + is the 

occupation fraction of state j at temperature T (300 K 
throughout this work), and r is an intrinsic relaxation 
time. The latter is a critical parameter, which we fix 
to for -1 = 1.6 meV, a value compatible with estimates 
obtained from measured DC mobilities [29j [30j [55] . The 
self-consistent potential is obtained from the external po- 
tential (/> ext and the induced charges as 
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FIG. 4: Definition of n in armchair and zigzag nano- 
triangles. We define n as the number of hexagons along the 
side (n = 2 in these examples). 



Appendix C: Size dependence 



where vu> is the Coulomb interaction between electrons 
at orbitals I and V . Finally, we calculate extinction cross- 
sections and near-field intensities from the self-consistent 
charge p = [x°/(l — v ■ x°)} ' ext derived from Eqs. 
(??) upon illumination by an external field E ext (i.e., 
we take 0® xt = —Hi • E ext , where Hi is the position 
vector of site I). For example, the cross section reads 
a abs _ (4 7ra; / c )l ni { Q /} 5 where a = (1/E ext ) J2 t x\p\ is the 
nanoisland polarizability. Further computational details 
can be found elsewhere l24l. 



We show in Fig. [5] extinction spectra for armchair nan- 
otriangles of two different sizes. The main effect of reduc- 
ing the size of the graphene island is an increase in the 
plasmon energies, which is significantly larger than the 
1/y/a dependence on side length a predicted by a clas- 
sical electromagnetism description [56]. However, like in 
classical theory, the absorption cross-section produced by 
these plasmons is roughly proportional to the area. 



2. Classical electromagnetic modeling 

We compare our quantum-mechanical results with 
classical finite-element (COMSOL) calculations in which 
the graphene is described as a thin film of thickness 
t and dielectric function 1 + Airier /out, where a(uo) = 
(ie 2 Ep /k% 2 ) /(uj + zr _1 ) is the Drude optical conductiv- 
ity. 

Appendix B: Geometry of the nanoislands 

Most of the calculations reported in this paper refer to 
graphene nanotriangles of either armchair or zigzag edges 
centered around the middle point of a carbon hexagon. 
We consider the edge atoms to be passivated with hydro- 
gen, and therefore, we treat them as the rest of the atoms 
in the structure, with the same hopping to their nearest 
neighbors. The number of atoms in the nanotriangles is 
3n(n + 1) for armchair edges and (n + 2) 2 — 3 for zigzag 
edges, where n is the number of hexagons along the side. 
For example, the n = 2 triangles of Fig. [4] contain 18 
and 13 carbon atoms, respectively. The side length of 
the triangles is a = (3n — 1) ao for armchair edges and 
a = \^3nao for zigzag edges, where ao = 0.142 nm is the 
nearest-neighbor distance. 



Appendix D: Collective character of the plasmons 

In the random-phase approximation (RPA) here used 
to describe the optical response of graphene nanoislands, 
the induced density is writen as p = [x° / (1— v-x°)]'^ ext in 
terms of the external potential </> ext , the Coulomb interac- 
tion between atomic sites v, and the noninteracting RPA 
susceptibility [24j (see Appendix [A]). The collective char- 
acter of the plasmons is clearly captured in a multiple- 
scattering fashion by the use of the self-consistent poten- 
tial [1/(1 — v - x°)} ' ext m the above expression. We 
compare in Fig. [5] these full RPA calculations (Fig. [5]3,d) 
with those obtained by removing the self-consistency and 
calculating the induced charge as p = \° -</> ext (Fig.pkf). 
This leads to dramatic changes in the plasmon spectrum: 
the plasmon dispersion with varying doping charge is re- 
moved, and the spectral features are now peaked at the 
positions of dominant electron-hole-pair transitions (see 
Fig. [3J3). This is in contrast to the plasmons obtained 
from self-consistent calculations, the energies of which 
depend on the interaction among different electron-hole- 
pair virtual transitions, as a manifestation of the collec- 
tive electron motion involved in these excitations. 
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FIG. 5: Size dependence and collective character of nanotriangle plasmons. We study plasmons in a graphene 
armchair nanotriangle of side a — 3nm (left), compared with those of the a — 7nm island considered in the this paper and 
reproduced here for convenience (right). We show the carbon atomic structure of the nanotriangles in (a) and (b), respectively. 
The absorption cross-section normalized to the graphene area is calculated for these two nanotriangles from the self-consistent 
RPA ((c) and (d)), compared with the non-self-consistent RPA ((e) and (f)). The spectral dependence of the cross section is 
given as a function of the number of elementary charges in the nanoislands, Q/e. 



Appendix E: Dependence on the choice of triangle 
center 

Plasmons in graphene nanoislands are very sensitive 
to small structural details of the carbon lattice. Besides 



the strong dependence on the type of edges observed in 
Fig. [2j we show in Fig. [6] a large variation in the plas- 
monic spectrum with the position of the triangle cen- 
ter. In Fig. [6^L,d,g and in the rest of the figures of this 
work, we present calculations obtained with the center 
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FIG. 6: Dependence on the choice of triangle center. Atomic structure (a-c), states near the Dirac point (d-f), and 
extinction spectra (g-i) of armchair nanotriangles with a — 7nm side length centered at the center of a carbon hexagon (H- 
centered) (a,d,g), at a carbon atom (C-centered) (b,e,h), or at the center of a nearest-neighbor carbon bond (C-C-centered) 
(c,f,i). 



located at the middle point of a carbon hexagon (H- 
centered). When the triangle center is chosed at a carbon 
site (C-centered, Fig.|6]3,e,h) or at the center of a nearest- 
neighbors bond (C-C-centered, Fig.[6j3,f,i), the corners of 
the atomic structure present a certain degree of asym- 
metry (see Fig. [6}d,c) and the one-electron energies are 
substantially modified (Fig.[6^,f). This leads to dramatic 
changes in the extinction spectra (Fig. [6ja,i) , which are 
still exhibiting a clear plasmon dispersion with varying 
doping charge. In particular, the lowest-energy plasmon 
of the C-C-centered triangle (Fig.[6j) shows a bump in en- 
ergy at zero doping, in contrast to the lowest-energy plas- 
mon of the H-centered triangle (Fig. [6^), which steadily 
increases in energy with doping. 



Appendix F: Effect of zigzag edges and zero-energy 
electronic states 



Zigzag edges are known to produce electron states near 
the Dirac point that are concentrated near the atomic 
edges. The number of such states is equal to the dif- 
ference between the number of atoms in the two carbon 
sublattices of the graphene nanostructure [33l |35j [37] . 
This is the case of the C-centered triangle of Fig. [6}d, 
which contains zigzag edges near the corners that lead to 
the existence of one zero-energy state. This state is half- 
filled in the neutral structure. Interestingly, adding or 
removing a single electron to the nanotriangle does not 
produce any significant change in the plasmonic spectrum 
because it ends up in an edge-localized electronic state 
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FIG. 7: Emergence of plasmons upon filling all zero-energy electronic states in a zigzag triangle. The zigzag 
triangle represented in (a) has 24 zero-energy states per electron spin at the Dirac point (see (b)), half of which are occupied 
in the neutral configuration. No significant change in the plasmon energy and strength is observed by adding electrons (holes) 
until the zero-energy states are fully occupied (empty). This happens for \Q\/e = 24. In contrast, the addition (removal) of 
each further electron beyond this point produces sizable changes, as shown in the extinction cross- section plotted in (c) for 
|Q|/e > 24. 



that cannot undergo large dipole transitions. However, 
the addition or removal of two electrons causes a sizeable 
shift. One can therefore speculate that electrons added 
(removed) to (from) the zero-energy states have no ef- 
fect in the collective plasmon excitations. Another piece 
of evidence pointing to this direction is provided by the 
zigzag triangle considered in Fig. [3]3,d,f, which contains 
a large number of zero-energy states and does not suffer 
any plasmon shifts with the addition of a small number 
of electrons or holes. In order to test this hypothesis, 
we consider the same zigzag triangle in Fig. [7| subject 
now to much higher doping. Interestingly, the plasmonic 
spectrum remains nearly unchanged while electrons or 
holes are added to the nanoisland. However, when the 
zero-energy states are completely filled (depleted) and 
further electrons (holes) are added to the structure, the 
plasmons undergo strong variations in energy with ev- 
ery new electron (hole). This is strong evidence that 



zero-energy states do not contribute to plasmons, and 
the tunability observed in graphene islands with vary- 
ing doping is quenched by the presence of these states. 
In summary, zigzag edges produce a net difference in the 
number of atoms in the two carbon sublattice of graphene 
nanoislands, leading to the presence of zero-energy states 
localized near the edges, the optical excitation of which 
is rather inefficient because they cannot undergo large 
dipole transitions, and therefore, zigzag edges are detri- 
mental to the tunability of graphene plasmons via doping. 



Appendix G: Effect of disorder in the edges 

The robustness of the spectra of armchair nanotrian- 
gles is tested in Fig. [8] against the degree of disorder in 
the edges. In particular, this figure shows that the spec- 
tral features of an immaculate doped a = 3nm trian- 
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FIG. 8: Effect of disorder in the edges. Spectral changes in doped armchair nanotriangles when disorder is introduced in 
the edges. The spectral features of a perfect nanotriangles of side length a = 3 nm (a) are increasingly split as additional atoms 
are added to the edges: 2 atoms in (b) and 4 atoms in (c). 



gle are increasingly split as extra atoms are added to 
the edges. Additionally, low-energy features are becom- 
ing more intense. Incidentally, the number of atoms in 
both carbon sublattices remains the same in all the struc- 
tures considered in Fig. [8] (i.e., the two atoms added in 
each step belong to different sublattices), so that there 
are neither zero-energy states nor effects associated with 
dummy zero-energy electrons, as described in the previ- 
ous section. 



Appendix H: Formation of the Dirac cone in large 
islands 

We show in Fig. [9] the formation of the Dirac cone 
structure as the size of the nanotriangle is increased. The 
Dirac cone emerges in the momentum-energy representa- 
tion of the electronic bands in extended graphene. How- 
ever, the momentum is not a good quantum number in 
finite nanoislands due to the lack of translational symme- 
try. Instead, we have taken the spatial Fourier transform 
of each and all of the one-electron states and represented 



the sum of their intensities as a function of wave vec- 
tor. This produces an energy-momentum structure that 
is reminiscent of the Dirac cone, which is more clearly 
emerging as the size of the nanotriangles increases. 



Appendix I: Maximum absorption cross-section 

The absorption cross-section at the peak of the fea- 
tures associated with the excitation of plasmons can ex- 
ceed the geometrical area of the graphene. We show this 
in Fig. [lO] for the plasmons of the island considered in 
Fig. [2^. These cross-sections are thus large enough to 
produce complete optical absorption in a surface deco- 
rated with nanoislands [39]. Importantly, the maximum 
cross-section is directly proportional to the intrinsic de- 
coherence time r, which can be increased due to ther- 
mal and disorder effects [28]. Nonetheless, our choice of 
for -1 = 1.6 meV is compatible with estimates based upon 
reported DC mobilities for high-quality graphene [28] . 



[1] S. Nie and S. R. Emory, Science 275, 1102 (1997). 

[2] K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, 

I. Itzkan, R. R. Dasari, and M. S. Feld, Phys. Rev. Lett. 

78, 1667 (1997). 
[3] P. Johansson, H. Xu, and M. Kail, Phys. Rev. B 72, 

035427 (2005). 

[4] C. E. Talley, J. B. Jackson, C. Oubre, N. K. Grady, C. W. 



Hollars, S. M. Lane, T. R. Huser, P. Nordlander, and 
N. J. Halas, Nano Lett. 5, 1569 (2005). 

[5] L. Rodriguez- Lorenzo, R. A. Alvarez-Puebla, 
I. Pastoriza- Santos, S. Mazzucco, O. Stephan, M. Ko- 
dak, L. M. Liz-Marzan, and F. J. Garcia de Abajo, J. 
Am. Chem. Soc. 131, 4616 (2009). 

[6] X. Qian, X.-H. Peng, D. O. Ansari, Q. Yin-Goen, G. Z. 



11 



probability density 



max 



> 1 
>^ 

§ -i 

-2 

-3. 





\ / 






rx 




: 




/ \ 

4 \ 







K 

a = 7nm 7V=816 
ac edges H-centered 



r r 




r 



a=14nm 7V=3168 
ac edges H-centered 



a = 28nm TV = 12870 
ac edges H-centered 



FIG. 9: Building up the Dirac cone structure. Intensity of the Fourier transform of all one-electron wave functions in 
armchair nanotriangles of increasing size. This is the continuation of Fig. [3^l towards larger sizes of the nanoisland. 




FIG. 10: Absorption cross-sections exceeding the 
graphene geometrical area. Peak absorption cross- 
sections associated with the plasmons of Fig. [5J:. 



Chen, D. M. Shin, L. Yang, A. N. Young, M. D. Wang, 

and S. Nie, Nat. Biotech. 26, 83 (2008). 
[7] L. Hirsch, R. Stafford, J. Bankson, S. Sershen, B. Rivera, 

R. Price, J. Hazle, N. Halas, and J. West, Proc. Natl. 

Acad. Sci. 100, 13549 (2003). 
[8] R. Asahi, T. Morikawa, T. Ohwaki, K. Aoki, and 

Y. Taga, Science 293, 269 (2001). 
[9] P. V. Kamat, J. Phys. Chem. B 106, 7729 (2002). 
[10] M. Grzelczak, J. Perez- Juste, P. Mulvaney, , and L. M. 

Liz-Marzan, Chem. Soc. Rev. 37, 1783 (2008). 
[11] R. Chapman and P. Mulvaney, Chem. Phys. Lett. 349, 

358 (2001). 

[12] T. Hirakawa and P. V. Kamat, Langmuir 20, 5645 (2004). 
[13] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009). 



[14] L. Ju, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, 
H. A. Bechtel, X. Liang, A. Zettl, Y. R. Shen, et al., Nat. 
Nanotech. 6, 630 (2011). 

[15] Z. Fei, G. O. Andreev, W. Bao, L. M. Zhang, 
A. S. McLeod, C. Wang, M. K. Stewart, Z. Zhao, 
G. Dominguez, M. Thiemens, et al., Nano Lett. 11, 4701 

(2011) . 

[16] H. Yan, X. Li, B. Chandra, G. Tulevski, Y. Wu, M. Fre- 
itag, W. Zhu, P. Avouris, and F. Xia, Nat. Nanotech. 7, 
330 (2012). 

[17] H. Yan, T. Low, W. Zhu, Y. Wu, M. Freitag, 
X. Li, F. Guinea, P. Avouris, and F. Xia, arXiv 0, 
arXiv:1209.1984vl (2012). 

[18] J. Chen, M. Badioli, P. Alonso-Gonzalez, S. Thongrat- 
tanasiri, F. Huth, J. Osmond, M. Spasenovic, A. Cen- 
teno, A. Pesquera, P. Godignon, et al., Nature 487, 77 

(2012) . 

[19] Z. Fei, A. S. Rodin, G. O. Andreev, W. Bao, 
A. S. McLeod, M. Wagner, L. M. Zhang, Z. Zhao, 
M. Thiemens, G. Dominguez, et al., Nature 487, 82 
(2012). 

[20] N. K. Emani, T.-F. Chung, X. Ni, A. V. Kildishev, Y. P. 

Chen, and A. Boltasseva, Nano Lett. 12, 5202 (2012). 
[21] Z. Fang, Z. Liu, Y. Wang, P. M. Ajayan, P. Nordlander, 

and N. J. Halas, Nano Lett. 12, 3808 (2012). 
[22] Z. Fang, Z. Liu, Y. Wang, P. M. Ajayan, P. Nordlander, 

and N. J. Halas, ACS Nano 6, 10222 (2012). 
[23] D. Pines and P. Nozieres, The Theory of Quantum Liq- 
uids (W. A. Benjamin, Inc., New York, 1966). 
[24] S. Thongrattanasiri, A. Manjavacas, and F. J. Garcia de 

Abajo, ACS Nano 6, 1766 (2012). 
[25] P. R. Wallace, Phys. Rev. 71, 622 (1947). 
[26] B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New J. 

Phys. 8, 318 (2006). 
[27] L. Brey and H. A. Fertig, Phys. Rev. B 75, 125434 (2007). 
[28] M. Jablan, H. Buljan, and M. Soljacic, Phys. Rev. B 80, 

245435 (2009). 

[29] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 
Firsov, Science 306, 666 (2004). 



12 



[30] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature 438, 197 (2005). 
[31] N. Armbrust, J. Giidde, P. Jakob, and U. Hofer, Phys. 

Rev. Lett. 108, 056801 (2012). 
[32] N. M. R. Peres, R. M. Ribeiro, and A. H. Castro Neto, 

Phys. Rev. Lett. 105, 055501 (2010). 
[33] J. Fernandez-Rossier and J. J. Palacios, Phys. Rev. Lett. 

99, 177204 (2007). 
[34] M. Ezawa, Phys. Rev. B 76, 245415 (2007). 
[35] A. R. Akhmerov and C. W. J. Beenakker, Phys. Rev. B 

77, 085423 (2008). 
[36] Z. Z. Zhang, K. Chang, and F. M. Peeters, Phys. Rev. B 

77, 235411 (2008). 
[37] M. Wimmer, A. R. Akhmerov, and F. Guinea, Phys. Rev. 

B 82, 045409 (2010). 
[38] M. Zarenia, A. Chaves, G. A. Farias, and F. M. Peeters, 

Phys. Rev. B 84, 245403 (2011). 
[39] S. Thongrattanasiri, F. H. L. Koppens, and F. J. Garcia 

de Abajo, Phys. Rev. Lett. 108, 047401 (2012). 
[40] Y. L. Luo, Y. S. Shiao, and Y. F. Huang, Nat. Nanotech. 

6, 247 (2011). 

[41] D. Subramaniam, F. Libisch, Y. Li, C. Pauly, 
V. Geringer, R. Reiter, T. MashofT, M. Liebmann, 
J. Burgdorfer, C. Busse, et al., Phys. Rev. Lett. 108, 
046801 (2012). 

[42] S. Kim, S. W. Hwang, M.-K. Kim, D. Y. Shin, D. H. 
Shin, C. O. Kim, S. B. Yang, J. H. Park, E. Hwang, 
S.-H. Choi, et al., ACS Nano 6, 8203 (2012). 

[43] Y. Kobayashi, K. ichi Fukui, T. Enoki, K. Kusakabe, and 
Y. Kaburagi, Phys. Rev. B 71, 193406 (2005). 



[44] J. Tian, H. Cao, W. Wu, Q. Yu, and Y. P. Chen, Nano 

Lett. 11, 3663 (2011). 
[45] Q. 6. Girit, J. C. Meyer, R. Erni, M. D. Rossell, 

C. Kisielowski, L. Yang, C.-H. Park, M. F. Crommie, 

M. L. Cohen, S. G. Louie, et al., Science 27, 1705 (2009). 
[46] H. Shu, X. Chen, X. Tao, and F. Ding, ACS Nano 6, 

3243 (2012). 

[47] C. K. Gan and D. J. Srolovitz, Phys. Rev. B 81, 125445 
(2010). 

[48] J. Wu, W. Pisula, and K. Mullen, Chem. Rev. 107, 718 
(2007). 

[49] X. Feng, M. Liu, W. Pisula, M. Takase, J. Li, and 

K. Mullen, Adv. Mater. 20, 2684 (2008). 
[50] F. Schedin, A. K. Geim, S. V. Morozov, E. W. Hill, 

P. Blake, M. I. Katsnelson, and K. S. Novoselov, Nat. 

Mater. 6, 652 (2007). 
[51] P. Bharadwaj, A. Bouhelier, , and L. Novotny, Phys. Rev. 

Lett. 106, 226802 (2011). 
[52] J. Kundu, F. Le, P. Nordlander, and N. J. Halas, Chem. 

Phys. Lett. 452, 115 (2008). 
[53] F. H. L. Koppens, D. E. Chang, and F. J. Garcia de 

Abajo, Nano Lett. 11, 3370 (2011). 
[54] A. Bostwick, T. Ohta, T. Seyller, K. Horn, and E. Roten- 

berg, Nat. Phys. 3, 36 (2007). 
[55] K. I. Bolotin, K. J. Sikes, Z. Jiang, M. Klima, G. Fu- 

denberg, J. Hone, P. Kim, and H. L. Stormer, Sol. State 

Commun. 146, 351 (2008). 
[56] J. Christensen, A. Manjavacas, S. Thongrattanasiri, 

F. H. L. Koppens, and F. J. Garcia de Abajo, ACS Nano 

6, 431 (2012). 



